Is there a linear potential at short distances? 
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We argue that the lattice data of the static potential can be explained by perturbation theory up to energies 
of the order of 1 GeV once renormalons effects are taken into account. 



One expects on general grounds the matching 
coefficients of effective field theories of QCD to 
suffer from renormalon ambiguities. This means 
that, in principle, one can not compute these 
matching coefficients with infinity accuracy in 
terms of the (short distance physics) parameters 
of the underlying theory. Being more specific, 
let us consider the effective Lagrangian in some 
schematic form: 

oo 

C = E — c *°s > (!) 

s=0 

where O s stands for the operators in the effective 
theory and c s for the short distance matching co- 
efficients. If we pick up one matching coefficient 
c(v/m), it would have the following perturbative 
expansion in a s : 

oo 

C (,/m)=c + ^^ +1 , (2) 
n=0 

where a s (in the MS scheme) is normalized at the 
scale v. Its Borel transform would be 

n=0 

and c is written in terms of its Borel transform as 

oo 

c = a + /d, e -'/". B H (i ). (4) 



The ambiguity in the matching coefficient reflects 
in poles 1 in the Borel transform. If we take the 
one closest to the origin, 

6B[c](t)~— (5) 
a — t 

1 In general, this pole becomes a branch point singularity 
but this does not affect the argumentation. 



where a is a positive number, it sets up the max- 
imal accuracy with which one can obtain the 
matching coefficients from a perturbative calcu- 
lation, which is (roughly) of the order of 

Sc~r n «af, (6) 

where n* ~ — . Moreover, the fact that a is posi- 
tive means that, even after Borel resummation, c 
suffers from a non-perturbative ambiguity of or- 
der 

al3 

5c~(A QCD )—. (7) 

In order to relate the matching coefficient with 
the observable, we use the effective field theory 
as a tool to provide the power counting rules in 
the calculation. Since we are relating one observ- 
able (a renormalon free object) with a matching 
coefficient suffering of some renormalon ambigu- 
ity, there must be another source of renormalon 
ambiguity as to cancel this one. The latter comes 
from the calculation of the matrix element in the 
effective theory. If the matrix element in the ef- 
fective theory is a non-perturbative object, the 
ambiguity of the matching coefficient is of the 
same size than this non-perturbative object that 
we do not know how to calculate anyhow. Nev- 
ertheless, if the matrix element in the effective 
theory is a perturbative object, the only way it 
has to show the renormalon is in a bad perturba- 
tive behavior in the expansion parameters of the 
effective theory (what is happening is that the 
coefficients that multiply the expansion parame- 
ters are not of 0(1) due to the renormalon, i.e. 
the renormalon is breaking down the assumption 
of naturalness implicit in any effective theory). 
This means that in the latter situation the ob- 
servable is less sensitive to long distance than the 



matching coefficient itself. Thus, if, for instance, 
we wanted later on to get the short distance pa- 
rameters from that observable, we are not doing 
an optimal job, since we use an intermediate pa- 
rameter (the matching coefficient) that can not 
be obtained with better accuracy than the ones 
displayed in Eqs. (6)-(7), whereas the observ- 
able is less sensitive to long distances (the same 
problem also appears if we want to relate two 
weakly-sensitive-to- long-distance-physics observ- 
ables through a more sensitive to long distance 
sensitive matching coefficient). 

The renormalon of the matching coefficients 
can be spurious (it is not related to a real non- 
perturbative contribution in the observable) or 
real (it is related to a real non-perturbative con- 
tribution in the observable). In fact, this distinc- 
tion depends on the observable we are considering 
rather than on the renormalon of the matching 
coefficient itself. The point we want to stress is 
that at the matching calculation level it makes no 
sense this distinction. Therefore, there is no ne- 
cessity to keep the renormalon ambiguity in the 
matching coefficients. Thus, our proposal is that 
one should figure out a matching scheme where 
the renormalon ambiguity is subtracted from the 
matching coefficients. This is the program we will 
pursue here for the specific case of effective field 
theories with heavy quarks and in particular for 
the singlet static potential. 

The static potential is the object more accu- 
rately studied by (quenched) lattice simulations. 
This is due to its relevance in order to under- 
stand the dynamics of QCD. On the one hand, 
it is a necessary ingredient in a Schrodinger-like 
formulation of the Heavy Quarkonium. On the 
other hand, a linear growing behavior at long 
distance is signalled as a proof of confinement. 
Moreover, throughout the last years, lattice sim- 
ulations [1,2] have improved their predictions at 
short distances allowing, for the first time ever, 
the comparison between perturbation theory and 
lattice simulations. Therefore, the static poten- 
tial provides a unique place where to test lattice 
and/or perturbation theory (depending on the 
view of the reader), as well as an ideal place where 
to study the interplay between perturbative and 
non-perturbative physics. This is even more so 



since the accuracy of the perturbative prediction 
of the static potential has also improved signifi- 
cantly recently [3-6]. 

Let us first review the status of the art nowa- 
days. The prediction for the perturbative static 
potential at two loops [3] indicated the fail- 
ure (non-convergence) of perturbation theory at 
amazingly short distances. This failure of per- 
turbation theory is not solved by the inclusion 
of the leading logs at three loops computed in 
[4,5] nor by performing a renormalization group 
(RG) improvement of the static potential at next- 
to-next-to-leading log (NNLL) [6]. On the other 
hand, it was soon realized that the static poten- 
tial suffered of renormalons [7] and that the lead- 
ing one (of O(Aqcd)) cancelled with the leading 
renormalon of twice the pole mass [8]. Never- 
theless, the prediction of perturbation theory for 
the slope of the static potential should not suffer, 
in principle, from this renormalon and could be 
compared with lattice simulations, which, indeed, 
only predict the potential up to a constant. This 
comparison was performed in Ref. [9] , where they 
compared the RG improved predictions (without 
including ultrasoft log resummation) of pertur- 
bation theory up to next-to-next-to-leading order 
with lattice simulations. They indeed found that 
the discrepancies with lattice simulations and the 
lack of convergence of the perturbative series re- 
mained. They also found that the difference be- 
tween perturbation theory and lattice could be 
parameterized by a linear potential in a certain 
range. Therefore, it seemed to support the claims 
of some groups [10,11] of the possible existence of 
a linear potential at short distances. Such claims 
contradict the predictions of the operator prod- 
uct expansion (OPE), which state that the lead- 
ing non-perturbative corrections are quadratic in 
distance at short distances 2 . 

On the other hand, it has been argued [12] 
that, within a renormalon based picture, pertur- 
bation theory agrees with the phenomcnological 
potentials aimed to describe Heavy Quarkonium. 
Moreover, In Ref. [13] (see also [12]), it has been 

2 This contradiction is indeed so only if one believes that 
perturbation theory can be applied for the shortest dis- 
tances available in lattice simulations ~ 1 — 4 GeV. We 
thank V.I. Zakharov for stressing this point to us. 



shown that perturbation theory can indeed re- 
produce the slope of the static potential given 
by lattice simulations at short distances by us- 
ing the force instead of the potential as the basic 
tool without the need to talk about renormalons. 

In these proceedings, we would like to review 
the results of Rcf. [14], which tried to clarify 
further the above issues and support the renor- 
malon dominance picture by comparing the po- 
tential computed in lattice with the RS static po- 
tential (a renormalon free definition of the poten- 
tial obtained in Ref. [15]). 

The energy of two static sources in a singlet 
configuration reads (in the on-shcll (OS) scheme) 



E(r) 



2m s + lim ~hi(W n ) 



(8) 



In the situation Aqcd *C 1/r, it can be computed 
order by order in a s (v) (y ~ 1/r) and in the 
multipole expansion (see [16,6]) 

E{r) = 2m s + ^,os(r,^us) + SE vs (r, i>us) • (9) 

mos is the pole mass. The static potential reads 
(a s = a s (v)) 



+i 



(10) 



n=0 



where V n = V n (r, v, v\js)- The first three co- 
efficients Vb, V\ and V2 are known [3]. The 
log-dependence on v of V3 can also be obtained 
by using the ^-independence of V^os- The log- 
dependence on j^us of V3 is also known [4] . There- 
fore, the only unknown piece of V3 is a f/^us - 
independent constant. Its size has been estimated 
in Ref. [15] assuming renormalon dominance. For 
rif — 0, it reads 



V 3 (r, 1/r, 1/r) = 1/r x (-76.1075) . 



(11) 



We will use Eq. (11) in the following for our 
estimates of V3. The static potential is v- 
independent. By setting v = 1/r, we could effec- 
tively resum the Invr terms. The RG-improved 
expressions would read 



V s ,os = J2v n (r,l/r,vus)a s (l/r) 



n+l 



(12) 



and we would have expressions for n = 0, 1, 2, 3. 
In the above expression we have considered 
lnr^us ~ 1- Since lnr^us 1 for some range 
of the parameters, one could be in the situation 
where one also has to resum ultrasoft logs as it has 
been done in Ref. [6]. Nevertheless, explicit cal- 
culations show that, at least for the range of pa- 
rameters studied in this paper, higher order ultra- 
soft logs are subleading even if sizable. However, 
for definitcness, we will work with the resummed 
expression (the physical picture does not change 
anyhow), which we will add to Vs(r, 1/r, 1/r): 



V 3 (r,l/r,iAjs) = V 3 (r,l/r,l/r) 



(13) 



+ T^ a l{r ^log — ) f 

6po \a s {v us ) 

The RS scheme was defined in Ref. [15] aim- 
ing to eliminate the renormalons of the matching 
coefficients appearing in Heavy Quarkonium cal- 
culations. In particular, the leading renormalon 
of the heavy quark mass and the static potential. 
We refer to Ref. [15] for details. Here, we just 
write the relevant formulas needed for our analy- 
sis. 

Analogously to the OS scheme, the energy of 
two static sources in a singlet configuration reads 



E s {r) = 2m RS (v f ) 



(14) 



+ ( lirn^ - \n(W n ) - 2dm RS (vf) 



where 



5m RS (v f ) = Y 5 ™RS, n a™ +1 

71 = 1 

^ T(n + l + b-k) 

x h Ck r ( 1+& - fc ) ' 



(15) 



n=0 



and rriRs = mos ~ <5wrs- If the beta function 
were known to infinity order in perturbation the- 
ory, it would be possible to obtain all the coeffi- 
cients b and Ck [18]. In practice, only b and co,i,2 
are known (see [18,19,15]). For N m only an ap- 
proximate calculation is possible along the line of 
Ref. [20]. This computation has been done in 



Ref. [15] (see also [21]). The result reads (for 
nj = 0) 

N m = 0.424413+0.174732+0.0228905 = 0.622036, 

where each term corresponds to a different power 
in u (in the Borel plane) of the calculation (see 
Ref. [15] for details). We see a nice convergence. 
This number will be the one we will use in the fol- 
lowing. This number should be equal to —2Ny, 
where Ny is the normalization factor of the renor- 
malon of the static potential. Ny was also ap- 
proximately computed in Ref. [15]. For rif = 0, 
it reads 



-1.33333+0.499433-0.338437 = -1.17234. 



In this case the convergence is not so good but 
we have an alternating series. In any case, we see 
that both values appear to be quite close: 



2N n 



Ny 



2N„ 



Ny 



= 0.059. 



(16) 



We take this as an approximate indication of the 
error in the evaluation of N m . 

By using Vrs, we expect the bad perturbative 
behavior that one finds with Vbs due to the renor- 
malon to disappear. Let us see that this is indeed 
so. We consider the OS and RS potentials at 
different orders in perturbation theory. We will 
work with v = 1/0.15399 r^ 1 . We display our 
results in Fig. 1. We see that we do not have 
the gap between different orders in the pertur- 
bative expansion (or it is dramatically reduced) 
in the RS scheme compared with the OS scheme. 
Although not displayed in these plots, if we had 
worked with v = 1/r, the RS results also converge 
to the same value. 

We can now compare with lattice simulations 
up to a constant in Fig. 2. We see that our results 
are convergent to the same potential, which, as we 
can see, corresponds to the lattice potential. 

We now study possible non-perturbative effects 
in the static potential. Any non-perturbative ef- 
fects should be small and compatible with zero 
since perturbation theory is able to explain lattice 
data within errors. We can make this statement 
more quantitative by using the lattice data ob- 
tained in Ref. [2] where the continuum limit has 
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Figure 1. Plot of toVrs(?") and roVos{r) °^ ^ ree 
(dashed line), one-loop (dash-dotted line), two- 
loops ( dotted line ) and three loops ( estimate ) plus 
the leading single ultrasoft log (solid line). For 
the scale of a s (y), we set v = l/0.15399r ~ 1 . 
v us = 2.5 r^ 1 and Uf = 2.5 r^ 1 . 



been reached. For these lattice points the system- 
atic and statistic errors are very small (smaller 
than the size of the points). Therefore, the main 
sources of uncertainty of our (perturbative) eval- 
uation come from the uncertainty in the value of 
A^g (±0.48 r^ 1 ) obtained from the lattice [17] 
and from the uncertainty in higher orders in per- 
turbation theory. We show our results in Fig. 3. 3 
The inner band reflects the uncertainty in Aj^<t 
whereas the outer band is meant to estimate the 
uncertainty due to higher orders in perturbation 
theory. We estimate the error due to perturba- 
tion theory by the difference between the NNLO 
and NNNLO evaluation. This is different from 
the procedure followed in Ref. [14], where we 
allowed Co, the three-loop coefficient of the static 
potential in momentum space (for the specific def- 
inition see [14]) to vary by by ±146. We now 
believe that this may overestimate the size of the 
perturbative errors, as a change of this magnitude 
should also be correlated with the change of the 
other variables of the computation (in particular 



3 We have performed the estimation of the errors with the 
results of Fig. 2 iy = constant), since they do not in- 
troduce errors due to the evaluation of the renormalon. 
Nevertheless, a similar conclusion had been achieved if the 
analysis had been done with v = 1/r. 




Figure 2. Plot of r {V RS (r) -V KS {r') + E latt .(r')) 
versus r at tree (dashed line), one-loop (dash- 
dotted line), two-loops (dotted line) and three 
loops (estimate) plus the leading single ultrasoft 
log (solid line) compared with the lattice sim- 
ulations [2] Ei a tt.{r). For the scale of a s (y), 
we set v = 1/0.15399 r^ 1 . v us — 2.5 r^ 1 and 
r' = 0.15399 r . 



Figure 3. Plot of ro(V R s(r) - V RS (r') + E latt .(r')) 
versus r at three loops ( estimate ) plus the leading 
single ultrasoft log (dashed line) compared with 
the lattice simulations [2] Ei att (r). For the scale 
of a s {v), we set v = 1/0.15399 r^ 1 . v us = 2.5 r^ 1 
and r' — 0.15399 r . The inner and outer band 
are meant to estimate the errors in A-^g and per- 
turbative. For further details see the main text. 



Nv s )- Moreover, it would not be consistent with 
the convergent pattern the series shows so far. In 
any case, the difference is not large as one can see 
by comparing Fig. 3 with Fig. 4, for which the 
errors has been obtained with this last method. 
The difference in the error band, not in the central 
value, of Fig. 4 with Fig. 9 and, accordingly, Fig. 
10 of Ref. [14] is due to the fact that we use here 
the 1-loop a s running for the variation, whereas 
in Ref. [14], we used the three- loop a s running 
(the use of the three-loop a s running introduces 
an incomplete set of subleading logs, therefore we 
prefer to stick ourselves to the strict results from 
the renormalization group method even if the in- 
troduction of these subleading logs would signif- 
icantly diminish the errors), plus a missprint in 
the formula used for the plot. 

The usual confining potential, SV = err, goes 
with an slope a = 0.21GeV 2 . In lattice units we 
take: a — 1.35 r^ 2 . The introduction of a lin- 
ear potential at short distances with such slope 
is not consistent with lattice simulations as we 
can see from Fig. 5. This is even so after the 
errors considered in Fig. 3 have been included. 



This should not come as a surprise since this lin- 
ear potential appears as an effect of long distance. 
Therefore, it follows that the use of the Cornell 
potential (with the perturbative static potential 
emanated from QCD instead of a pure 1/r po- 
tential times a fitted constant) as a phenomeno- 
logical fit of the static potential introduces sys- 
tematic errors if the typical inverse Bohr radius 
scale of the heavy quarkonium system to study 
lye in the short distance regime as it is the case, 
for instance, for the T(IS'). 

On the other hand, recently, there have been 
claims about the possible existence of a linear po- 
tential at short distances [10,11,22]. Expected it 
to be of different physical origin than the long 
distance linear potential, it may have a differ- 
ent slope than the (long distance) static poten- 
tial discussed above. It would be very important 
to discriminate its existence, since such behav- 
ior at short distances is at odds with the OPE. 
This is indeed possible in some cases, since the 
short-distance linear potential expected in Ref. 
[10,11] have an slope of the order of magnitude 
of the long-distance confining potential that we 




Figure 4. Plot of r (V RS (r) -V RS (r') + E latt .{r')) 
versus r at three loops ( estimate ) plus the leading 
single ultrasoft log (dashed line) compared with 
the lattice simulations [2] Ei att (r). For the scale 
of a s {v), we set v = 1/0.15399 r§ . v us = 2.5 r^ 1 
and r' = 0.15399 r . The inner and outer band 
are meant to estimate the errors in A^jg and cq. 
For further details see the main text. 



Figure 5. Plot of ro(V RS (r)-V RS {r') + Ei att .(r') + 
o~(r — r')) versus r at three loops (estimate) plus 
the leading single ultrasoft log compared with the 
lattice simulations [2] Ei att (r). For the scale of 
ct s {v), we set v = 1/0.15399 r^ 1 . v us = 2.5r ~ 1 
and r' = 0.15399 r . The dashed band is meant 
to estimate the combined error due A^g and per- 
turbation theory (see Fig. 3). 



have already ruled out above. Therefore, we can 
conclude that no linear potential exists at short 
distance (with the present estimates for its slope) . 
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